Solution of the Unanimity Rule on exponential, uniform and scalefree networks: 
A simple model for biodiversity collapse in foodwebs 
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We solve the Unanimity Rule on networks with exponential, uniform and scalefree degree dis- 
tributions. In particular we arrive at equations relating the asymptotic number of nodes in one of 
two states to the initial fraction of nodes in this state. The solutions for exponential and uniform 
networks are exact, the approximation for the scalefree case is in perfect agreement with simulation 
results. We use these solutions to provide a theoretical understanding for experimental data on 
biodiversity loss in foodwebs, which is available for the three network types discussed. The model 
allows in principle to estimate the critical value of species that have to be removed from the system 
to induce its complete collapse. 
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I. INTRODUCTION 



Unanimity rule (UR) is generally associated with mod- 
els, where a (binary) state of an atom or agent can change 
only if all of its direct neighbors are in the other state, re- 
spectively. Usually UR is formulated in a network frame- 
work, where a node becomes 'activated' only if all the 
nodes pointing to it - through directed links - are 'acti- 
vated'. These models have attracted some recent interest 
because of a number of important real world applica- 
tions. Maybe the most relevant example is the modeling 
of biodiversity based on foodwebs. Species are nodes in 
this web. If one species is food for another species this 
is indicated by a directed link in the foodweb, pointing 
from the eaten species to the eater. Usually one species 
does not depend on a single other species but in gen- 
eral has a more diversified menue. In the picture of the 
foodweb this means that each node is pointed at from 
several other neighbors. Imagine that all species in a 
hypothetical foodweb exist in one of two states: alive 
or extinct. If all the neighbors pointing to a particu- 
lar node i, are extinct the node itself has no more food 
to live on and will go extinct in the next timestep as 
well; the relation to the UR becomes apparent. UR has 
been studied experimentally in actual foodwebs to model 
biodiversity loss [U [2] ■ Earlier efforts on modeling food- 
web topology and its fragility have been conducted with 
respect to clustering, fragmentation, robustness, and de- 
gree distribution under the assumptions of random and 
non-random extinction of species [3] . In [4] small world 
effects are investigated and it was concluded that food- 
webs are not random networks with Poissonian degree 
distribution. Niche models - as a method of sampling 
surrogate foodwebs - are studied in |6j [7] which (in 
the low connectancy limit) display robust scaling prop- 
erties of foodwebs which is in good agreement with data 
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from several field studies, see e.g. |5]. Further, their con- 
clusion that the degree distribution of foodwebs decays 
exponentially rather than scale-free is in good agreement 
with [TJ [5]. Studies concerned with the robustness of 
foodwebs generally use the UR as an update mechanism, 
propagating initial extinction of species. 



Surprisingly, the dynamics of Unanimity Rule, which 
is a generalization of the Majority Rule of opinion dy- 
namics [SI EH EH HH H3] and reminds on features of the 
Voter model [H E03 US E] , the Axelrod model [HI US] 
as well as of Boolean networks [201 El] , is poorly known 
[2"2l |2"5] and has been put on a more mathematical ba- 
sis only recently |24j . The unanimity model as presented 
there can also be viewed as a limiting case of a model 
for decision making scenarios [25]. Note that UR dif- 
fers from these previous models by the fact that it is 
irreversible, i.e. once a node has reached the activated 
state, it remains in it. This irreversibility of UR makes it 
an excellent candidate not only for modeling biodiversity 
as above but also for the adoption of new technologies, 
such as MMS [26 , by interacting customers. Techno- 
logical standards are generally irreversible once they are 
adopted by a population. Another specificity of UR is the 
fact that it is purely deterministic, i.e. once the topology 
of the underlying network is fixed and an initial number 
of nodes are activated, the entire dynamics is determined. 
In contrast, the Voter model, when applied to complex 
networks, involves a random step when a node chooses 
an interaction partner among its neighbors. Similarly, in 
the Majority Rule, a node choses randomly between two 
nodes among its neighbors to form a majority triplet. 



In this paper we start from the dynamical equation for 
UR on networks and solve it analytically for the special 
cases of exponcncntial and scale free networks. We then 
compare the results with experimental findings obtained 
from actual foodwebs [T] [5] . 
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FIG. 1: First two steps of UR starting from an initial network of 7 nodes, 2 of them being activated. Initially there is only 
one node among the non-activated nodes that satisfies the unanimity rule. It gets therefore activated at the first time step. At 
that time, there is a new node whose 2 incoming links come from activated nodes. It gets activated at the second time step. 
This system gets fully activated at the fourth time step. In the context of foodwebs for example node 2 eats species 1 and 4, 
and serves as food for species 3 and 5. A black node means that it is extinct. (From [24]). 



II. UNANIMITY RULE 

UR is implemented amongst agents in a network com- 
posed of N nodes connected through L directed links. 
Each node exists in one of two states: activated or inac- 
tivated. The number of nodes with indegree i (number 
of links pointing to it) is denoted by Ni and depends on 
the underlying network structure, i.e. the indegree dis- 
tribution is given by pi — Ni/N. It is a fixed quantity 
that does not evolve with time. The fraction of activated 
nodes at time t we denote by at , for the fraction of acti- 
vated nodes with an indegree i, we write a i t - Initially (at 
t = 0) there is a fraction of do activated nodes. Obviously 

a t = '^2pia l ,t (1) 

i 

holds. The unanimity rule is defined as shown in Fig. 1. 
If all the links arriving to an inactivated node i originate 
from nodes which are activated at t— 1, i gets activated at 
t. Otherwise, it remains inactivated. At each time step 
every node is considered for an update, the dynamics is 
synchronous. The process is iterated until the system 
reaches a stationary, frozen state, characterized by an 
asymptotic value afi na i = doc . The UR problem is to pre- 
dict (Zoo from the knowledge of the do and the structure 
(indegree distribution) of the network. 

To solve problems of this nature a two step approach 
was suggested [2H first map the problem onto an 
update equation, second find the asymptotic solutions 
of the latter. To derive the update equation, note that 
the probability that i randomly chosen nodes are initially 
activated, is a l (i is an exponent). The average number 
of nodes with indegree i and which respect the unanimity 
rule is therefore Noq. Among these nodes, Niaoa l were 
already activated initially. This is due to the fact that 
the total number of nodes with (in)degree i which are 
initially activated is NidQ. Consequently, the number of 



nodes that gets activated at the first time step is 

A i)0 = (Ni - N iao ) 4 , (2) 

and, on average, the total number of occupied nodes with 
indegree i evolves as 

Ai,i = A i<0 + A ifi ■ (3) 

At the next time step, the average number of nodes with 
indegree i, which respect the unanimity rule and which 
are outside the initial set is (Ni — Nao)a\. Among those 
nodes, A^o have already been activated during the first 
time step, so that the average number of nodes which get 
activated at the second time step is 

A itl = (Ni-N iao )(ai-ai) . (4) 

Note that Eq. Q is valid because no node in A^i also 
belongs to A i0 - This is due to the fact that each node 
can only be activated by one combination of i nodes in 
our model, so that no overlap is possible between A^i 
and Aifl. By iterating it is straightforward to show that 
the contributions A i t read 

A M = (JV i -JV i a )(t4-oj_ 1 ) , (5) 

with a_i = 0, by convention. The number of activated 
nodes evolve as Ai^t+i — Ai, t + A^t, and by dividing by 
Ni, one gets the equations for the fraction of activated 
nodes a, S [0,1] 

<k,t+i = OM+(l-ao)(oi-aJ_x) = a;,o+(l-aoK ■ ( 6 ) 

Now a t is the convex sum of Eq. ^ with weights accord- 
ing to the indegree distribution pi, i.e. a t = ^2iPiO,i,t- 
Finally we get for the update equation 

at+i = a + (1 - a ) ^Pia\ . (7) 

i 

Numerical solutions for t — ► oo are found in Fig. [2] for 
several different indegree distributions pi (symbols). 
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III. EXPONENTIAL AND SCALEFREE 

Let us focus on the special choices of the exponential 
(pi oc e~ Xt ), and scalefree (pi oc i~ x ), indegree distribu- 
tions, to analytically uncover the behavior of Eq. ([7]). 
The simpler cases, pi = b~i\ and pi — 8 i2 have been solved 
in 

For the exponential case we have to find an expres- 
sion for the term X^iP^t; where M is the maxi- 
mum indegree, i.e. M < N. With the exact identity, 



c(l — c )/(! — c), the exact asymptotic equa- 



tion for the exponential case is found to be 



ao + (1 - a )a o 



1 



1 - e" AJM 1 - e- A aoo / ' 

where we first have set c = exp(— A)aoo for the non- 
normalized contribution of Pi<i\, and second c = exp(— A), 
to account for the norm of the distribution pi. Finally 
we have taken the limit t — > 00. 

We get two important limiting cases: First, the uni- 
form distribution is recovered as the limit A — > 0, 



a + (1 - a ) 



1 



M \ 1 



,M 



(9) 
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FIG. 2: Asymptotic solutions, a x , of the UR on three types 
of networks vs. initial condition ao. Symbols represent the 
numerical solutions of Eq. |7f. The networks used for this 
solutions have been generated with A = — 1, M = 10 for the 
exponential, M = 4 for the uniform and A = —0.5, M = 10 
for the scalefree case. Lines are the theoretical results of Eqs. 
(8}, (9| and <[TTJ , for the same parameters. The critical value 
is where the plateau at 1 is reached (all species extinct). 



where we have used the fact 
exp(-MA))/(l - exp(A)) = M. 
Second, the large system limit M 



that liniA^o(l — 



00 is 
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ao + (1 - a )a c 



1 



,-A 



1 — e~ A a„ 



(10) 



The scalefree case is treated similarly. We approximate 

J2tLi *~ A a* « 1 + lt/2 1/2 dx x ~ Xa t and E^i i~ X ~ 

1 + Jy^ 1 ^ 2 dx x~ x . Defining the two-sided incomplete 

gamma-function as j(a,b,x) = f^dy y x ~ 1 exp(— y) the 
asymptotic solution for the scalefree unanimity rule reads 



floo = a Q — (l—a )X- 



^oo| A 7( 



|ooo|,(Af +|)|a 0O |, -A) 



[M ■ 



r A -(|)- A -A 



(ii) 

where A = A — 1. The quality of the results for the 
exponential, uniform and scalefree cases is seen in Fig. [2j 
where the asymptotic value is plotted against the initial 
ao. Points represent the numerical solution to Eq. ([7]) for 
realizations of networks of the characteristics specified in 
the figure caption. Solid lines are the analytical results 
from Eqs. ([I]), ^ and (111. These equations show that 
the larger the exponents A (for exponential and scalefree) 
the larger the critical value (where the plateau at 1 is first 
reached) becomes. Similarly, the larger M the larger the 
critical value. 



In [TJ |5] experimental results on foodwebs are pre- 
sented, which are the outcomes from numerous field stud- 
ies, see References 9,11,15,28-37 in [2j. In p] the number 
of species having died out as a result of the initial removal 
of a given fraction of the total population is presented. 
This constitutes a UR on networks, where - in our nota- 
tion - the fraction of activated links a t are the ones that 
have died out after t timesteps. The reported networks 
there have been exponential, uniform and scalefree. We 
have re-drawn the data from [1] for two examples for each 
,network type [27] in Fig. [3] (symbols). Unfortunately, the 
network parameters can not get reconstructed from [2J. 
The data presented contains in- and outdegrees mixed 
together, which makes it impossible to estimate A and 
M for the indegree only (which is needed for the UR). 
The lines represent our theoretical results where we have 
estimated the parameters in the following way. For the 
exponential we fixed A to the values reconstructed from 
Fig. 2 in [2J, and varied M until good fit was obtained. 
Of course M opt is less than M reported in Fig. 2 of [2\. 
The uniform case turned out to be optimal with M opt 
about half of the reported value, the scalefree case was 
fitted after fixing M to about half of the reported case 
and varying A. 
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FIG. 3: Experimental data from [HE] (and references therein), for the number of asymptotically died out species on measured 
exponential (a), uniform (b) and scalefree (c) distributions (for randomly removed species). Lines are curves from our theoretical 
solutions with parameters as explained in the text. 



V. DISCUSSION 

We have solved the UR on exponential, uniform and 
scalefree networks. These solutions relate the structure 
of the underlying network together with an initial loss of 
diversity, clq with the final diversity in the system. Once 
the structure of the underlying network is known, these 
solutions allow further to find the critical value of initial 
'species removal' a§ rlt at which (and beyond which) a 
total collapse of the system will appear. The dynamics on 



foodwebs (who eats whom) is a particular example of this 
UR, the reported underlying networks are of exponential, 
uniform and scalefree nature. We have shown that not 
only do our results compare well with simulations but 
also with actual experimental data. The message of this 
work is to point out the potential danger of uncontrolled 
anthropogenic species removal. 
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